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Abstract. Red horizontal-branch (RHB) stars have been selected from the Hipparcos catalogue to investigate 
their kinematics and spatial distribution. Hipparcos parallaxes, literature radial velocities and Hipparcos proper 
motions, together with models for the gravitational potential of the Milky Way allow a calculation of the actual 
velocity vectors and the orbits of the RHB stars. The velocity characteristics are used to define a halo population 
sample (HPS) in the collection of RHBs. The orbits lead statistically to an overall z-distance probability distri- 
bution, showing that the RHBs exhibit two populations, a disk one having a scale height of ftdisk — 0.6 kpc and a 
halo one of ~ 4 kpc. We have investigated the influence on our results of parallax accuracy and of a demarcation 
line in the HRD between the RHB and the red-giant (RG) star region. Neither of them show marked effects. We 
have performed the orbit analysis using the potential model of Allen & Santillan as well as of Dehnen & Binney. 
The results differ only slightly for the disk population, showing that these potential models are not a critical part 
of such orbit investigations. RHB scale height values are smaller than those found earlier for sdB stars, most likely 
because the samples of stars used had different spatial distributions a priori. The data do not allow us to specify 
a trend in the kinematic behaviour of star types along the horizontal branch. 
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1. Introduction 

While it is rather simple to study the structure of other 
stellar systems from images alone, the true shape of the 
Milky Way is very hard to determine. With the help of ever 
growing galactic surveys one is finally able to make use 
of highly accurate positions and motions of thousands of 
stars. In conjunction with mass models of the Milky Way, 
such a wealth of data gives us the first detailed glimpses 
of the structure and evolution of our own galaxy. 

Choosing a single stellar type, in particular one with 
well-defined luminosity, allows us to work with clean sam- 
ples. Late main-sequence (MS) stars were used to probe 
the solar vicinity within several hundred parsecs, while 
giants and A type MS stars were observed up to two kilo- 
parsecs away. For regions of negligible interstellar extinc- 
tion the s tellar density field was der ived from the se sur- 
veys (e.g. lBeckejll94(ll lElviuslll95lL and llJpgrenl ITflBai . 
Results from these studies established the idea of a multi- 
component galaxy, consisti ng of a central bulge , a thin 
disk, a thick disk (found bv lCilmore fc ReiJl983(l and an 
encompassing halo. 

Further information on the exact shapes of these com- 
ponents can then be used to constrain theories of galactic 
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evolution. Models for the formation history of the disk are 
separated into two families. If the thin disk is formed af- 
ter the thick disk, one speaks of a "top down" scenario. 
Different exam p les of thi s mod el type are descr i bed i n 
iGilmore ('l984"), 'Sandage ("199^), iNorris fc Rvanl l)l99l|) . 
and Burkert et al., (1992') . If the thin disk is formed be- 
fore the thick disk, it is called a "botto m up" s c enario . 
Members of this g r oup are the models of Norris fl987L 
' on et all l)l993|) . iKroupa (,200 2) and Gilmor e fc 
(19 83). A summ ary of the various models is given by 
lMaiewskill)l993^ . 

With the astrometry satellite Hipparcos, an auto- 
mated survey of stellar positions with unprecedented ac- 
curacy was made possible. Some 120000 stars were ob- 
served for the main mission. Data were collected in the 
Hipparcos Catalogue (esa 199^. In addition, the satel- 
lite provided two-colour photometry for slightly more than 
one million stars. The data were reported in the Tycho 
Catalogue (|ESAIIi 997|) . 

Among the brighter old stars easily identifiable from 
optical observations are the RR Lyr stars and blue sub- 
dwarf stars. These are stars of the horizontal branch (HB), 
consisting of a ~ 0.5 M0 He core surrounded by a hydro- 
gen shell ranging in mass from negligible to up to about 
0.5 M0. From thin to thick shells, the stars are classified 
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Fig. 1. CMD of single HiPPARCOS stars with A{B - V) < 
0.025 and — < 20%. The stars inside the window are 
our selected RHB candidates. The window borders are: 
{B-V)> 0.5, My e {-0.5 1.8}, Mv < -2.8 + 7-(B-F), 
My < 7.3 - 6.2 • {B - V). The dash-dotted hne marks the 
window separation border of Sect. 14. 4.!^ while the dashed 
line marks the alternative cut-off of Scct. l4.4.'2l 



as sdO, sdB, HBB and HBA stars, followed by the RR Lyr 
stars, and finally the "massive" , red HB (RHB) stars. Of 
these, the sdB and RR Lyr stars are the easiest to spot 
and characterize. The former show a simple relation be- 
tween the line strengths of the Balmer series and Tcff and 
logg. The latter, as variable stars, are identified by their 
characteristic light curve. 

Several studies exi st deal i ng wit h the spatial distri- 
bution of sdB stars. iHebeJ l)l986l) and iMoehler et alJ 
found vertical distributions with scale heights 
200 PC, way too small for old stars. However , 



1199! 

of 



IVilleneuve et all l)l995b|) and IVilleneuve et all l)l995a|) 

showed that larger values are more realistic . In two exten- 
sive studies of the kinema tics of sdB stars ( de Boer et all 
ll997llAltmann et al.l20ol it could be shown that the sdB 
stars seen in the solar neighbourhood are divided in two 
populations, a disk one with a scale height of ~ 1 kpc and 
an extended one chracterized by a scale height of ~ 7 kpc. 
Furthermore, trends were seen in the kinematics of HB 
star samples. The cooler Blue HB (BHB) star samples 
seemed to be more do minated by the halo popul ation than 
the sdB star sample ('Altma nn fc de Boeill2000l) . 

In this study we investigate the kinematics and spatial 
distribution of the red HB stars. Finding RHBs in a sky 
field is not without difficulty. By just looking at colours, 
one can easily confuse them with MS stars, subgiants, red- 
clump stars and blue-loop stars. Without knowledge of 
absolute magnitudes it is therefore not possible to select 
RHBs just by colour. Using Hipparcos data we are able 



Table 1. Numbers of stars at various selection steps as 
a function of relative parallax error. The three BB2k 
columns show: number of RHB sample stars found, total 
number of measurements, number of stars that are fainter 
than 7.3 mag, the completeness limit of the HiPPARCOS 
mission. 



Att 

TT 


HIP 


RHBs 


stars 


BB2k 
meas. 


>7.3 mag 


< 10% 


17295 


547 


455 


480 





< 15% 


28700 


1393 


828 


890 


28 


< 20% 


38114 


2241 


1044 


1127 


103 


< 25% 


45622 


3084 


1215 


1313 


211 


< 30% 


51231 


3795 


1327 


1437 


305 



to resolve this matter. Moreover, it gives us a rather large 
sample of RHBs to work with. 

To characterize the large scale distribution of stars of 
a given type one in principle would need to know the 
location of such stars in complete large distance sam- 
ples. For our HiPPARCOS RHB stars this is illusionary, 
since the sample is complete out to perhaps 100 pc (see 
Sect . 13.2.11) . However, the kinematics of the RHB stars also 
give information about past and fut ure locations. We will 
calculate the galacti c orbits (as in Ide Boer et al.l Il997t 
lAltmann et al.ll2004l) of the RHB stars and analyze the 
orbit statistics to arrive at mean spatial distributions. 

The organisation of this paper is the following. We 
begin with the assembly of our data sample in Sect.|21 
Errors and selection effects of the sample are described in 
Sect.Ol The analysis of the orbits of the selected RHBs 
is given in Sect. 01 while Sect. El compares results using 
different models for the galactic potential. We conclude 
with a final discussion in Sect.El 

2. Establishing the RHB data sample 

iTautvaisiend iimd) and iRosd l|l985|) give observational 
methods to find RHB stars through photometriy and spec- 
troscopy. However, the simplest way to define the state of 
evolution of stars is given by their location in the HRD. 

Using HiPPARCOS data we created a colour mag- 
nitude diagram of stars not listed in the Catalogue 
of the Components of Double and Multiple Stars 
toomma nget fc NvslEoO^ . having {B - V)< 0.025 mag, 
and relative parallax error smaller than 20%. The RHB 
stars could then be selected from the diagram by choos- 
ing an appropriate range in My and (B — V) by eye 
(see Fig.nj . CMDs of globular clusters like 47 Tuc helped 
choose the proper magnitude and colour ranges. We also 
confirmed that the spectroscopically identified RHB stars 
from Behr (2003,) mostly fulfil the selection criteria. 

For more on the complexities of the boundary line to- 
ward the giant branch see Sect. 13. 2."^ 

Five samples were assembled from stars within this 
window. These sets differ in parallax precisions, with each 
having ^ lower than 10%, 15%, 20%, 25%, and 30%, 
respectively. Thus, with such a criterion, we almost get a 
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volume-complete sample for spheres with a rather blurry 
boundary. Effects from this sample definition are discussed 
in Sect. OH for the 10%, 20%, and 30% samples. TableHl 
gives the number of stars in the samples at different points 
of the assembly process. 

To reconstruct the full spatial motion, radial veloci- 
ties are required. HiPPARCOS did not have a dedicated 
instrument to measure Urad and so we chose the recent 
rel ease of the General Catalog of mean radial velocities 
bv lBarbier-Brossat fc FigonI (|2D0dl); (hereafter BB2k). By 
cross-checking our candidate RHBs with BB2k we found 
radial velocities for 35% up to 83% of our HIP stars, de- 
pending on the parallax accuracy (Table^. This already 
hints at selection effects that are discussed in Sect . 13. 21 



3. Error limits and selection effects 

3.1. Error limits 

Measurement errors for the Hipparcos objects vary not 
only with brightness but also with celestial position. This 
is due to the scanning law because certain areas in the 
sky were observed more often than others. Therefore, for 
all HIP stars of Hp < 9 mag (the faintest of our stars 
has Hp — 8.86 mag) we only cite median precisions. 
They are 0.97 mas for the parallaxes and 0.88 mas/yr 
and 0.74 mas/yr for /i* = /j,q, • cos S and fig, respectively. 
Systematic errors were estimated to be below 0.1 mas. 

Errors from BB2k depend on the source from which the 
radial velocity was taken. Single high resolution spectra, of 
course, yield better results compared to big surveys at low 
resolution. However, errors are usually below 10 km s~^ 
and never higher than 20 km for our samples. 

Hipparcos proper motion uncertainties for stars of 
our sample are equal to ~ 2 km s~^. Therefore the total 
error in space velocity is dominated by the uncertainty of 
the radial component. Thus, we estimate the total error 
in space velocity to be about 10 km s~^. 

3.2. Possible selection effects 

3.2.1. Hipparcos stars and measurement accuracy 

Hipparcos was an astrometric mission that used a list of 
objects t o be observed , the H ipparcos Input Catalogue 
(HIC) bv lTuron et aTl (|l992h . A magnitude-complete set 
down to 7.3 mag was assembled, some 52000 objects. Any 
stars fainter than that limit came from individual propos- 
als. The selection effects so introduced cannot be quanti- 
fied. For the number of stars in our sample fainter than 
7.3 mag see Tabled 

Hipparcos parallaxes are usually determined at an 
average precision of 0.8 mas for high ecliptic latitudes, 
and around 1.2 mas near the ecliptic. High accuracy ar- 
eas for ^* are roughly at \/3\ > 40° (0.8 mas yr~^, down 
from 1.3 mas yr^^). For on average the standard error 
is 0.8 mas yr" "*^. No position-dependent colour errors are 
given in .ESAl l|l997l) . 



The mission's Hp passband was more sensitive to blue 
light (Fig 1.3.1 in IeSAI Il997|) . Redder stars, therefore, 
show a higher uncertainty in their measurements. Since 
our stars are usually brighter than V = 8 mag, this effect 
is small and probably only shows up in the crowded region 
of the RGB cut-off. 



3.2.2. Boundary of RHB and Red Clump 

The definition of our RHB sample in the CMD is a selec- 
tion effect in its own right. While the range in luminosity 
and the boundary on the side of the MS can be easily de- 
fined, the separation from the red-giant branch (RGB) is 
less certain. 

In particular, that region of the CMD contains the 
"Red Clump" , a region where three main types of stars 
are found. First, there are the RG stars evolving upwards 
in the CMD (but having an uneven distribution along the 
RG). These are mostly contributing to the red part of 
the Red Clump, and are thus excluded from our sample. 
Second, there are the red HB stars distributed horizon- 
tally toward the RGB but of which the ones with mass 
up to 1.5 Mq are located above and some what bluer than 
the reddest end o f the HB range (see e.g. ISweigartlll987L 
ISeidel et al .119871 for the overall shape of the HB). Finally, 
the yet more massive blue loop stars may contribute to the 
more luminous part of the Red Clump re gion. A full dis - 
cussion of the complexities can be found in| lGallartllll998l) . 

Since all HB stars have essentially the same life time 
irrespective of their mass, the possible presence of younger 
(and thus more massive) RHBs cannot be neglected. 
Young RHBs exhibit disk kinematics and we may expect 
the disk population to be more dominant in our sample. 

Effects of the choice of the boundaries in the CMD (see 
Fig.lll) will be discussed in Sects. 0^ and l4.4.3l 

3.2.3. Radial velocity catalogues 

The BB2k data are a compendium of measured radial ve- 
locities from various sources published before December 
1990. By looking at how many of the RHBs compared to 
the whole HIP stars can be found in the BB2k catalogue, 
we can estimate how well, compared to other stellar types, 
our RHB stars are represented in BB2k. 

We limit the search to objects brighter than 7.3 mag in 
V for completeness reasons. A total of 15752 single stars fit 
this criteria, with 2052 falling in our RHB window. BB2k 
lists entries for 9446 (60%) and 1101 (54%) of these cate- 
gories, respectively. We conclude that our RHB candidate 
stars are not significantly undersampled in BB2k. 

3.2.4. Overall effects 

Considering the properties of the catalogues used and our 
selection criterion, we can estimate the volume complete- 
ness of our sample. Hipparcos measured all stars down 
to 7.3 mag. In conjunction with our lower absolute mag- 
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Fig. 2. The Toomre diagram of our ^< 30% sam- 
ple. The circles show the peculiar velocity ?;pcc = 
^J^"^ + + {Q - QlsrY with respect to the LSR. Up 
to a pecuhar velocity of 50 km s~^ the sample is dis- 
tributed rather evenly. Upwards of these values the stars 
fall short in orbital velocity. These objects are halo and 
thick disk stars. 

nitude limit of 1.8 mag gives (assuming negligable extinc- 
tion) a completeness radius of 126 pc. The BB2k is not 
magnitude complete. Some 1600 RHB candidates with 
^<30% and brighter than 7.3 mag are not hsted in 
BB2k. A merging of the two catalogues thus shrinks the 
completeness radius. For ■^<30% the new limiting mag- 
nitude is 6.4 mag correspond i ng to a distance of 83 pc. 
In contrast to lAltmann et al.l l(2004l) , we are not limited 
towards the bright magnitudes. Thereby it is certain that 
all nearby candidate RHBs are included in our sample. 

For all objects with V < 7.3 mag we know the HIC to 
be complete. Measurement errors of these objects are such 
that we can confirm the correct CMD position for 91.3% 
of our RHB sample. The only other factor of influence is 
the BB2k, which gives Vrad for ~54% of the selected HIP 
stars. A total percentage of completeness for these bright 
stars is then found to be 91.3%-54%~50%. For our fainter 
objects this value decreases to ~20% if we assume that 
the HIC only contains half of all fainter RHBs. 

4. The spatial distribution of RHB stars 

4.1. Distribution from position and velocity 

A histogram of the circular velocity Q of our sample shows 
a high narrow peak located at roughly 6lsr- These are the 
thin disk stars. After taking out this group, the presence 
of a broader distribution centered around O ~ 180 km s~^ 
becomes apparent. It may represent the group often called 
the thick disk. Finally, there are stars covering a wide 
range of values of 8 G [—100,-1-340] km s~^. These are 
the halo stars in our sample. However, we cannot assign 
an individual star to these groups and therefore no velocity 
dispersions were calculated. 

The separation of the thick and thin disk stars can be 
seen clearly in a Toomre diagram (Fig.[5J). A large number 
of objects cluster in the region with Vpec< 50 km s"^. 
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Fig. 3. The Bottlinger diagram of our ^< 30% sample. 
The panels at right show the enlarged regions marked by 
squares in the left panels. Heavier points are for the Halo 
Population Sample (HPS) stars. The circle symbols denote 
solar values. Note the two groups of stars in the more 
detailed 6, <& panel. 



These are mainly thin disk stars that are moving along 
with the LSR. Stars in a range of 50 <Wpec< 150 km s""'^ 
are mostly thick disk objects. They are lagging behind in 
O because they are kinematically heated up compared to 
thin disk stars. This effect is known as the asymmetric 
drift of the thick disk (see e.g. lBinnev fc Trem ainc IQSj). 

All other stars, upwards of Vpcc— 150 km s^^, must 
be halo objects. Some have large perpendicular velocities, 
both with prograde and retrograde orbits. Stars with large 
peculiar velocities include those staying in the disk but 
having highly non-circular orbits. An example of these is 
the star with 8 ~ 340 km s~^, but which never strays 
further than 530 pc from the galactic plane. We make a 
rough selection of halo objects for further analyses and 
will refer to these 25 stars as the Halo Population Sample 
(HPS). 

The separation of the populations can be further in- 
vestigated with the help of a Bottlinger diagram (Fig.O. 
Stars of the above defined HPS are marked by heavier 
points and show velocities that are very different from the 
LSR. Other stars outside of the enlarged frames (Fig.O 
right side) are most certainly thick disk stars. 

The right half of Fig.|31 shows the values similar to the 
LSR in more detail. The 8, ^-diagram shows a separation 
into two groups. Such structures are well known and have 
been analysed in various studies. Several stellar streams 
were already found i n the U, V distribution of Hipparcos 
data. ISkulian et alJ l)l999|) even associated close star clus- 
ters with some of these streams. We refer to'Skuli an et alJ 
((1999i) and .Nordstrom et al.. (,20041) for more on that topic. 
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Fig. 4. Mer idian plots of som e of th e calculated orbits in 
a model bv lAllen fc SantillanI l|l99lh . The top row shows 
orbits of HPS stars. The middle row shows the orbit of the 
fastrunner star HIP 98269 (left) while the other two are 
thick disk candidates. The bottom row depicts thin disk 
orbits. 



4.2. Population assignment from orbit data 

We have calculated the orbits of our stars based 
on th e galactic potential model of lAllen fc SantillanI 
lll99lh Thereafter AS) using an orbit code developed 
bv lOdenkirchen fc Bros"ch3 ill992>l. a nd further a dapted 
for our work jde'^oer^ alJ 1997, Alt mann et al ] 12004). 
Orbits are calculated in steps of 1 Myr over a time span 




Fig. 5. Logarithmic histograms, that show the complete 
^< 30% sample (gray) and the contribution of our HPS 
stars (black). We have added a constant offset for visi- 
bility reasons. The ecc- histogram shows a separation of 
populations around 0.45. As expected, the halo stars show 
much higher eccentricities than the disk stars. In the nze- 
histogram the halo stars do not only build the tail of the 
distibution, but also account for low values. The peak on 
the right side of the diagram includes the stars with higher 
normalized z-extent than plotted. 
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Fig. 6. Diagram of O (gray triangles) and Omcd (black 
squares) against total kinetic energy and median kinetic 
energy, respectively. Filled symbols are stars of the HPS. 
The parabolas are isovelocity lines for V}_ , which is orthog- 
onal to Q. Medianised values lie on an almost straight line 
pointing from the LSR data point towards lower i?kin,mcd 
and emod- 



of 10 Gyr (since we are not modeling the evolution of the 
Galaxy, this long span, although physically unrealistic, is 
allowed for a good sampling of the shape of wide orbits). 
A selected few orbits can be seen in Fig.0] Most of the 
orbits are of a regular box shape. A few are irregularly 
shaped, perhaps a result of the close approach of the star 
to the galactic center and thus the interaction with the 
Bulge potential. 

By looking at histograms of the orbit eccentricities and 
normalized ^-extents f^ Fig. | ^, we qualitatively come to the 
same conclusions as in I Altmaiin et all l)2004l) . In the ecc- 
histogram, a separation between a high and low eccentric- 
ity sample can be made with the help of the HPS. It is 
roughly at ecc = 0.45, above which the halo stars domi- 
nate the distribution. In the nze-histogram, the stars at 
nze > 0.3 are clearly of the HPS. But HPS stars also con- 
tribute to the peak at low normalized z-extents, since our 
definition of the HPS does not exclude highly eccentric, 
coplanar orbits. 

Fig.El shows the kinetic energy versus orbital velocity 
of our star sample. Both the measured and medianised 
values are given. Looking at the median values, one can 
see that most of the disk stars are situated along a line 
pointing from the LSR towards lower E'kin.mcd and 0med- 
A clear gap can be seen at 9 110 km s~^, below which 
there are only HPS stars. 

Using the information from Fig.|21we have separated 
the HPS stars from the disk stars, a separation confirmed 
through the kinematic representation of Figs.EJEl andEJ 
However, individual stars of one group may still be in- 
cluded in the other group. 
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Table 2. The fit results of the z-probabihty distribution 
from Fig.|3 Scale heights and number densities are given 
for each fitted component. For a two-component fit no 
"thick disk" values were determined. 



fit constants 


2 comp. 


3 comp. 


/ithin [kpc] 


0.1 


0.12 


iVthin(O) " 


5.5-10^ 


4.8-10'^ 


/ithick [kpc] 




0.58 


iVthick(O) " 




51000 


fthalo[fcpc] '' 


5 


4.2 


iVhala(O) 


1100 


900 



° Values of A'^(O) come from the numerics of the orbit statistics 
and cannot be used as number of stars per unit volume. 
Halo values are based on three stars only (see Sect . 14.31 . 



4.3. Spatial distribution from orbit statistics 

From the sum of the orbits we obtain the z-probability 
distribution of ou r sample stars, as introduced in 
Ide Boer et al.l l|l997j) and further refined in lAltmann et all 
l)2004|) . We analvsed it for disk and halo populations with 
a weighted x^"fit. The large range in N{z) made a weight- 
ing scheme necessary that could take into account both the 
disk region and halo region data points. We chose 1/N{z) 
as a weight to accomplish that goal. A simple exponential 
relation of the form N(z) = N{0) ■ e^^l^ was fitted for 
each component, where iV(0) is the star number density 
in the local galactic plane and h is the scale height. Both 
sides of the distribution were used for the calculations. 

Analyses of the orbits showed that the data points with 
z > 8 kpc are from three stars only. They are HD 126778, 
HD 175305 and HD 184266. All the other HPS stars show 
strongly elliptic orbits close to the galactic plane. 

Fig.El shows the z-probability distribution, as well as 
the best fits. First, a two-component fit, representing the 
disk and the halo, was made. But in the range of 2-5 kpc 
the full z-distribution is not very well approximated (see 
left half of Fig.d). We therefore also made a three compo- 
nent fit. The results give a much better fit to the data than 
with only two components (right half of Fig.[7J|. Assuming 
that exponential functions apropriately describe the z- 
distribution, our data suggests the existence of more than 
two spatially distinct stellar groups. 

The results of the fits are given in Table|21 We note 
that the scale heights of the halo are derived from only 
three stellar orbits and are only listed for completeness. 

The results obtained are based on what will be called 
the base-line model. Below we will explore star selections 
with different parallax accuracies or selected with differ- 
ent magnitude and colour limits (as in Fig.^. We will 
also explore (in Sect .13 the effects of different models for 
the mass distribution in the Milky Way on our base-line 
results. 
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Fig. 7. The z-probabihty distribution of our ^< 30% 
stars with a two-component and a three-component z- 
distribution exponential fit, shown in the left and right 
halves of the diagram, respectively. Fitted curves are: thin 
dashed ones the individual population components (see 
Sect . 14. 3(1 ■ thick dashed ones are the complete fits. Note 
that the three-component fit is a better approximation of 
the actual distribution. 

4.4. Effects of varying the chosen sample structure 

4.4.1. Samples with 10,20,30% parallax accuracy 

In Sect .13 we briefly gave some details on star samples with 
different ^ constraints. By analysing the 10% and 20% 
samples we can derive some information on the spatial 
distribution of the stellar populations. This implies dealing 
with smaller samples and thus a deterioration in number 
statistics. 

Diagrams derived from the current velocities are qual- 
itatively the same for all the samples. A closer inspection 
of the Toomre diagram shows that a majority of our ob- 
jects are stars with thin disk orbit characteristics and only 
a minor fraction is of a halo like type. Disk stars are in- 
cluded in all the samples in a large enough quantity for 
analyses. In the U, ^-diagram, the bimodality of the data 
can still be seen, even in the 10% sample. This is due to 
the fact that the separation is observed for disk stars, and 
they are adequately represented in all the percentage sets. 
In contrast, only three of the HPS stars have a relative 
parallax error below 10%. This low number implies that 
with these we cannot make a meaningful analyses for the 
halo. For — <20%, the number of HPS stars is twelve. 

The orbital data thus confirms the similarity of the 
composition of the error-selected samples. For the disk 
stars, the histograms are of similar shape and HPS stars 
of the 10% and 20% sample are also present with eccen- 
tricities larger than 0.4. We conclude that no kinematical 
preference can be found in the closer samples. 

Even though we have very low numbers of stars for a 
z-distribution analysis in the 10% sample (see Table^, 
we give rough estimates on the scale heights. They are 
0.12 kpc for the thin disk, 0.34 kpc for the thick disk and 
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8 kpc for the halo. Scale heights of the 20% sample are 
0.12 kpc, 0.46 kpc, and 9 kpc, respectively. These results 
are not different in an essential way from those of Tabled 

4.4.2. Influences of the RHB-RGB separation line 

Another, perhaps critical aspect of our RHB selection from 
HiPPARCOS data is the heavily populated location of the 
window border at the RG branch. By excluding more of 
the critical red stars, we can check for influences of giants 
and see if our above results still hold. 

A new separation was adopted, by shifting the line to- 
wards brighter My by 0.5 mag (see dashed line in Fig.QJ. 
We thereby excluded almost half of the full ■^<30% sam- 
ple ending up with some 1942 Hipparcos stars. From 
these stars, 731 had 789 measurements documented in 
BB2k (1327 stars in the fuU sample, see Tab.d. 

With this new and restricted sample we noted several 
things. The velocity histogram shows that some of the 
stars lagging behind the disk rotation are now missing. 
The histogram of thick disk velocity stars now consists 
of several peaks. Stars with thin disk velocities are also 
less numerous, but the peak is still very similar in shape. 
Looking at the Toomre diagram allows the identification 
of the "missing" stars. We found that only 15 of the HPS 
stars are still included in the restricted sample, but the 
relative number of HPS to disk stars stays the same (1.9% 
compared to 2.1% for the base-line sample). In addition, 
a reduction of star numbers in the range of 100 km s^^< 
v± <150 km s^^ was noticed. The absence of these objects 
explains new gaps in the ecc-histogram (at ecc = 0.35) 
and 8 , i?kin-diagram (at 8 = 175 km s~^). In contrast, 
the Bottlinger diagram and the histogram of normalized 
z-extent are qualitatively the same as with the base-line 
sample. 

Thus, we derived the scale heights for this restriced 
sample via a x^-fit from the z-probability distribution 
results for the disks in /ithin = 0.12 kpc and /ithick = 
0.46 kpc. For the halo we get /ihaio = 4.2 kpc. 

4.4.3. Non HB stars 

In the selection window there is a potential contamination 
of stars that do not belong to the HB, such as massive core 



helium burners or RG stars (see Sect. 

Core Helium burners are more massive than RHB stars 
and evolve faster. The probability to find these inside our 
window is small. 

RG stars of solar metallicity found within our selec- 
tion window are v ounger than 1.25 billion years (Yonsei- 
Yale isochrones bv lKim et alJl20o3) . Assuming a constant 
Star Formation Rate (SFR) and an observed RGB only 
of solar metallicity, we would expect that an RGB of 
[^]= —0.5 dex would lie roughly 0.25 mag to the blue. 
RG Stars of that metallicity develop faster by s ome 50% 
(deriv ed from Yonsei-Yale evolutionary tracks bv iYi et'al] 
l2003l) . so this metal-poor RGB would in the CMD be at 



two-thirds of the observed RGB's density. However, the 
density bluewards of the observed RGB is much lower than 
that. Had the SFR in the past been smaller, then the con- 
tamination of our sample is smaller, too. If the SFR had 
been larger in the past, we would see many more RGs 
of bluer colour than actually are seen in the Hipparcos 
CMD. All this shows that a contamination of our RHB 
sample (if present) comes from young stars and not from 
old, metal-poor RG stars. 

We check for this contamination of our RHB sam- 
ple by separating the selection window in two parts. The 
Blue Sample is defined as having [B — V)< 0.85 mag, 
while the other stars are assigned to the Red Sample (see 
Fig.^. Stars that belong to the bluer sample are certain 
to be RHBs. These 99 stars include 6 HPS stars. The Red 
Sample contains 1228 stars, including 19 HPS stars. 

The analyses of the two samples suggest that there are 
no significant differences in kinematic behaviour. The Blue 
Sample does have a comparatively higher fraction of halo 
stars, including the two stars that have orbits reaching 
out to 2; > 10 kpc. Scale heights of the Blue Sample are 
found to be /ithin,biuo = 0.1 kpc, /ithick.biue = 0.8 kpc, 
and /ihaio,biuc — 13 kpc. For the Red Sample we obtain 

/ithin,rcd = 0.12 kpc, /ithick.rod = 0.46 kpC, and /lhalo,rod = 

1.8 kpc. The latter value is so low because the far-out 
moving halo stars are part of the Blue Sample. 

Had our base-line sample been contaminated by stars 
of non-HB nature, in other words mostly younger stars, 
then the z-distribution from the various samples should 
show significant differences in the thin disk component. 
This was not observed and we therefore conclude that we 
are essentially working with just RHB stars. 

5. Effects of galactic mass models 

Several models have been proposed to describe the 
true gravitationa l pote nti al of our gal a xy, e. g., by 
AS. iFlvnn et all lll996tl . IJohnston etal] l|l995|) . and 
iDehnen fc BinnevI ( 19981) (hereafter DB) . Thusfar we have 
used the AS model, with wq = 8.5 kpc and 8 = 
220 km s^^. Here we want to compare our results for 
the orbit statistics with those derived using one of the 
DB models. They are using different shapes for the disk 
and spheroidal potentials. We traced the orbits of our 
selected RHB candidates in their potential called "2b". 
It was chosen, because with vjq — 8.5 kpc and 8 — 
231 km s^^ it is the one most similar to the AS model. 
Fits to the models yield the following Gort constants: 
A = 12.95 km s^^kpc"! B = -12.95 km s^^kpc-i for 
the model by AS and A = 13.8 kms^^kpc^^ B ~ 
— 13.3 km s~^ kpc~^ for the "2b" model by DB. To further 
facilitate comparisons with previous results, we chose the 
same time constraints (10 Gyr backwards in 1 Myr steps). 

In Fig.|Slwe provide the meridional plots of the same 
stars as in Fig.^jfrom the "2b" potential. The orbits of the 
first two objects shown appear to have changed consider- 
ably in the meridional plot, but they still retain almost 
the same eccentricities and normalized z-extents. For the 
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Fig. 8. Meridian plots o f some of the calculated orbits 
in the "2b" potential bv IPehnen fc BinnevI l|l99i^ . The 
same stars as in Fig.0] were used to allow comparisons. 
The orbits shown are the ones differing most between the 
AS and DB potential models (see Sect. 13). 

following two, reaching far out in vo, ecc and nze change. 
Since these stars stay closer to the galactic center in the 
"2b" potential of DB, the ecc decreases, while nze in- 
creases. Other orbits for disk stars appear to be largely 
unchanged on the plotted scales. 

Differences in ecc and nze between the orbits in the 
AS and DB potentials for the whole 30% ensemble are 
given in Fig. El On average, a slight decrease of ecc can be 
seen in Fig.|ni(a), while in normalised z-extent (Fig.|ni(b)) 
the disk stars lie slightly above zero, documenting larger 
values of nze in the "2b" potential. A closer look at the 
data revealed this to be a simultaneous effect of an increase 
of Zmax and decrease of ^(zmax)- These differences arise 
from the different characteristics of the two potentials. 

Medianised values of O and iJ^jn are found to be sim- 
ilar as in Fig.0 Again, strong clustering along a roughly 
straight line is seen, which reaches from the model's LSR 
rotational velocity of 231 km s^^ towards lower energies. 
The gap in that line around 8 = 115 km is also seen, 
but it is not as wide as for the AS data. 

The z-probability distribution of our RHB candidate 
stars in the "2b" potential of DB is given in Fig.llUI 
Differences to the fit parameter values for the AS potential 
only arise for the halo. Its lower scale height and higher 
number density are, however, in line with the change in 
the shape of the orbits (compare Fig.0]with Fig.lHJ. 

6. Comparing the distribution of RHB stars with 
other types 

The results for the spatial distribution of the RHB stars 
can be compared with the distributions found for other 




0.8 1.2 

nzeAs 

Fig. 9. Comparison of the eccentricities (a) and normal- 
ized z-extents (b) of our stars in the potential of Allen & 
Santillan (AS) and the "2b" potential of Dehnen & Binney 
(DB). The HPS stars are given as filled symbols, the oth- 
ers as grey open circles. Orbits are less eccentric in the DB 
than in the AS potential. In nze there is a large scatter in 
the differences. 



stellar types. I n particular a compar ison with the distri- 
bution of sdBs ijAltmann et al.ll2004t) is of relevance. 



2-distribution - 
single components - 
full three component fit - 
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Fig. 10. The z-distribution of our ^< 30% RHB sample 
in the "2b" potential by DB. The maximum \Z\ values are 
roughly 18 kpc, which is lower than in a potential from AS. 
This flatter distribution of the halo group is reflected in 
the fit parameters: /ithin — 0.12 kpc; iVthin(O) — 4.8 • 10^; 
/ithick = 0.58 kpc; 7Vthick(0) = 51000; /ihaio = 3.4 kpc; 

iVhalo(O) = 1700. 
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6.1. General comparison 

In Table[21we have collected the data from our fit to the 
^-probability distributions of the 30% sample. We give 
values for all three components, the thin disk, thick disk, 
and halo. While the iV(0) values cannot be compared to 
real midplane densities, their ratios can still be used for 
this discussion. 

For the thin disk, iMaiewskil lIlQO,*^ gives in his re- 
view a lower scale height limit of 325 pc, a v alue mostly 
from successful starcount model s. Recentl v . IChen et all 
l|200l|) found 330 ± 3 pc, while IDu et alJ l|2003|) found 
320±15 pc. Our value of 120 pc is muc h smaller than those. 
However, there are measure ments bvlKuii ken fc Gilmor3 
l(l989l) . lKrouDa et alJ l| l993l). and lOiha et all l|l994l) that 
give 249 pc, 250 pc, and 240 pc, respectively, as scale 
height for the thin disk. It is not obvious why all these 
results differ so much. Since we are using calculated orbits 
to determine the scale heights, instead of just positions, 
the potential mod el may also come into play. 

The review bv IMaiewskil il993') gives thick disk scale 
heights of 1.3 to 2.5 kp c. Some more rec entlv published 
values are 0.58-0 .75 kpc l|Chen et al.l200ll). 0.64±0 . 03 kp c 
i Du et allliool . and 0.8-1.2 kpc l|Kerber et all l200ll) . 
Thus, the scale height of the thick disk is found to be in 
the range of 0.6 to 2.5 kpc. Our value of 0.58 kpc for both 
potentials is at the lower end of the above range. Since the 
more recent studies tend to find values of 1 kpc or less, we 
think that we are in line with those publications. 

Midplane thick-to-thin-disk numbe r ratios in the lit- 
erature are given from as low as 2% (lGilmore|ll984|) up 
to 13% ijChen et al.ll200ll) . There is, due to the way such 
data are fitted, an anti-correlation between density ratio 
and scale height. Our ratio of about 1% is lower than the 
values given above with a rather low scale height. The 
ratio of halo to thin disk stars taken from our data is 
M,.in(0)/i ythin(0) - 0.02% (- 0. 04% for the "2b" po- 
tential bv lPehnen fc Binnevlll998j) . Ty pically, these val- 
ues range fro m 0.0 5% to 0.4% (see e.g. IChen et ani200lL 
iKerber et al.ll200il) . These discrepancies are probably due 
to the small volume in space, from which we drew our 
data, the small sphere around the Sun. Other studies are 
based on selected fields with, in general, more distant 
stars. These are less likely to include thin disk stars, while 
thick disk stars are less likely to be included in our sample. 

Could this aspect (samples of different distance ranges) 
explain all the differences? A comparison between a spher- 
ical and conical survey of stellar populations with the same 
densities and scale heights gives hints to the answer. For 
a sphere the number of disk stars in the volume increases 
with the observed distance. This increase is much greater 
than for a narrow cone that looks away from the disk. 
For example, with our fit results from Sect. 14.31 we obtain 
a ratio of A^haio/A'thin — 0.17% for a sphere with radius 
2 kpc, while the same ratio for the cone of the same length 
towards a galactic pole is 5.5%. For a distance of 5 kpc 
the numbers are 0.33% and 43%, respectively. This ex- 
treme case should emphasise, that even with our very low 



midplane density ratio, a conical survey can find a large 
number of halo stars. Effects arising from the shape and 
depth of the surveyed volume are therefore important. 

Another, more fundamental, reason for the different 
results comes from the method of analysis. While some 
other studies explicitly assign the stars to one of the three 
stellar populations (thin disk, thick disk, halo), we opted 
not to do so. Our results are not influenced by errors in 
such assignments. Rather, our derived scale heights and 
number densities are only a result of the fitting procedure 
and are dependent of the weights used. Differences coming 
from the fit procedure are thus not negligible. We chose 
inverse counts as weights, because the fit function approx- 
imated the distribution equally well for both the disk and 
halo component. 

6.2. Are there differences in population along the HB? 

In lde Boer et alJ ()l997() the z-probability distribution was 
introduced as a tool to determine scale heights of stellar 
populations. The la test results on BHB sta rs involving this 
method are given in lAltmann et~al] l|2004^ From some 120 
sub-dwarf stars, they derived a scale height of 0.9±0.1 kpc 
for the thick disk. Another population with scale height of 
about 7 kpc was found in the data. These values are larger 
than the ones derived from our RHB sample (0.58 kpc and 
4.2 kpc) by almost a factor of two. 

For thi s paper the same orb it calculation tools were 
used as in lAltmann et all l)2004() . Any systematic effects 
(too restrictive potentials as in Sect. lH.l|l arising from the 
analyses should therefore be the same for both studies. 
But the fitting of the 2-distributions was done differently. 
We used the whole range in Z to do a weighted fit for 
three components simultaneously. 

However, the values are taken from data of two posi- 
tionally very different samples. The RHB sample is located 
near the sun, while the sdBs are generally farther away at 
high galactic latitudes. Consequently, from the RHBs we 
find ^-distributions that lack a well defined halo, while for 
the sdBs a thin disk component was barely seen. 

The large difference for the halo scale heights from the 
RHBs results from this "undersampling" in our sample. 
Another effect of this is that we get the lower midplane 
halo-to-disk number density ratios mentioned in Sect. 16. ll 

A greater number of far-out halo stars could surely be 
found in the HiPPARCOS and BB2k data, but we deciced to 
limit the observed volume by our already rather large rel- 
ative parallax error of 30%. Gaia (Ferryman ct al. 20o3), 
ESA's future astrometric mission will measure three- 
dimensional positions and velocities down to 17.5 mag. 
With such a wealth of data, we will be able to derive 
firmer results on the halo scale height of RHBs. 

The thick disk components from both studies are well 
represented, showing a large number of objects. Scale 
height differences in this population can therefore not be 
explained by low number statistics. Again, the selection of 
the samples influences the results, albeit at a lower level 
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than for the halo component. Whether the scale heights 
can be explained by an intrinsic difference in kinematics 
of RHBs compared to sdBs, cannot be decided. For that 
we would need to observe the same volume of space. 

6.3. Summary of the results 

1) The position and velocity diagrams (Figs.[21and|31) show 
that the RHB sample is a mix of disk stars (with Q ~ 200 
km s~^) and stars with lower orbital speed and thus larger 
perpendicular velocity, establishing the group of the asym- 
metric drift. A relatively small subset exhibits velocity and 
kinetic energy of a nature typical for the halo. 

2) The orbit shapes indicate a predominance of rather 
circular orbits. Again, a small subset exhibits very elliptic 
orbits in part reaching to large distances from the disk. 

3) The orbit statistics make clear that only a few of the 
HiPPARCOS-BB2k stars reach far into the halo. This does 
not imply that stars found to have highly excentric orbits 
but staying in the disk are not part of the halo population. 

4) No significant dependence of the thin disk parameters 
on the volume surveyed was found. However, the thick disk 
and halo show varying parameters due to the low number 
of objects. 

5) Different models of the distribution of mass in the 
galaxy do not produce significant differences in the scale 
heights of our sample. 

6) The vertical scale height of our RHB star samples and 
the ones of the earlier investigated sdB stars are signifi- 
cantly different. This either points at a different history of 
the progenitors of these groups or, given the widely differ- 
ent spatial sampling, demonstrates that source selection 
leads to strong biases with concomittant differences in the 
resulting average parameters. 

7) The range of thick disk scale heights, as derived from 
orbits as well as other methods, is either due to sample 
selection problems or is inconsistent with a simple galactic 
structure. 
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